Proximity effects and characteristic lengths in ferromagnet-superconductor structures 
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We present an extensive theoretical investigation of the proximity effects that occur in Ferro- 
magnet/Superconductor (F/S) systems. We use a numerical method to solve self consistently the 
Bogoliubov-de Gennes equations in the continuum. We obtain the pair amplitude and the local 
density of states (DOS), and use these results to extract the relevant lengths characterizing the 
leakage of superconductivity into the magnet and to study spin splitting into the superconductor. 
These phenomena are investigated as a function of parameters such as temperature, magnet polar- 
ization, interfacial scattering, sample size and Fermi wavevector mismatch, all of which turn out to 
have important influence on the results. These comprehensive results should help characterize and 
analyze future data and are shown to be in agreement with existing experiments. 
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I. INTRODUCTION 



The importance of understanding the characteristic 
length scales and geometrical effects inherent to het- 
erostructures consisting of ferromagnets in electrical 
contact with superconductors has received a consider- 
able reinforcement from the ever increasing advances in 
nanofabrication technology. These advances have made 
it possible (see e.g. Refs. |l|-^) to fabricate high quality 
nanostructures involving ferromagnets, as well as normal 
metals, and superconductors. In parallel, there has been 
significant progress in the development and refinement 
in tunneling spectroscopy techniques. Scanning tunnel- 
ing microscopy (STM) allows one to locally probe the 
electronic density of states (DOS) of hybrid systems over 
atomic length scales with sub-meV sensitivity. 

When a normal metal is in good electrical contact with 
a superconductor, superconductivity is weakened in the 
superconductor and induced in the normal metal. When 
the normaLmetal is not magnetic, this phenomenon is the 
traditionalo proximity effect which is described quantita- 
tively via the pair amplitude, F(r), which encompasses 
the spatial dependence of pair correlations in the both 
the superconductor and the normal metal. If the non- 
superconductor is a ferromagnet, the superconducting 
proximity effect is drastically modified by the finite ex- 
change field. Furthermore, the magnetic material can 
induce spin polarization in the superconductor, resulting 
in a magnetic proximity effect. The study of the spatial 
variation of both the pair amplitude and the local DOS 
are fundamental to the understanding of these nanostruc- 
tures. 

When considering such inhomogeneous systems, there 
are multiple length scales involved that must be eluci- 
dated. At T = the phase coherence in a clean nor- 
mal metal in contact with a superconductor decays in- 
versely with distance from the interface, with axharacter- 
istic length, £jv(0), that is essentially infinite.cl At finite 
temperatures, the phase coherence decays exponentially 



over a much reduced distance £at(T)B Conversely, at low 
temperature, the pair correlations in the superconductor 
become depleted near the interface over a length scale 
given by the zero temperature superconducting coherence 
length £0, while for temperatures close to T c , Ginzburg 
Landau theoryH predicts that the depletion is governed 
by a length scale that diverges at T c . Although 

the essentials of the standard proximity effect have been 
well understood for a long time.pl the length scales in 
the intermediate temperature regimes have not been sys- 
tematically or consistently studied for a normal metal- 
superconductor bilayer system, although self consistent 
microscopic calculations exist forJayered structures J3 and 
results within the quasiclassicatl description have been 
obtained. Other predictions are limited by being based 
on phenomenological or non-self consistent approaches. 

If the nonmagnetic normal metal is replaced with a 
ferromagnet (F/S junctions and structures), the rele- 
vant length scales in the problem are altered significantly. 
Naively, one would expect that all phase coherence would 
be lost in the magnet, since the superconductor and fer- 
romagnet have opposite types of long-range ordering: a 
ferromagnet favors parallel spin alignment and acts as 
an effective pair-breaker, while a superconductor is com- 
prised of Cooper pairs with (in the ordinary s-wave pair- 
ing considered here) antiparallel spin alignment. How- 
ever, a stable superconducting state can arise in the fer- 
romagnet in which ±hc Cooper pairs have a net center 
of mass momentum.EI The spin splitting in the ferromag- 
net introduces a new length scale £2 set by the differ- 
ence in the spin up and spin down Fermi wavevectors, 
£2 oc (fcp-f — kpi)^ 1 , which is typically much smaller than 

An interesting manifestation of this effect is the n 
phase junction comprised of a ferromagnetic material 
sandwiched between superconductors.E3~E3 This partic- 
ular interplay of ferromagnetism, arid .superconductivity 
has been studied for some timeJlJilj The peculiar os- 
cillatory state (originally introduced in the context of a 
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new superconducting .state that arises in when magnetic 
impurities are presentESE-j) leads, in the sandwich geom- 
etry, to a nonmonotonic dependence of the critical tem- 
perature on the ferromagnet layer thickness. E£rO Other 
worksc3 have focused on the variation of the Josephson 
current with temperature. For certain values of the ex- 
change field, spontaneous currentsc3 may arise in F/S 
heterostructures. 

For heterostructure configurations in which any of the 
material thicknesses are of order of or smaller than the 
largest of the relevant intrinsic lengths, size effects will 
play a role, and the finite thicknesses of the layers be- 
come important geometric lengths in the proximity ef- 
fects. It is then clearly preferable to tackle the prob- 
lem using a theory which does not involve coarse grain- 
ing over atomic length scales. It has been shown that 
in thin superconductor-normal metal bilayers, the inter- 
layer resistance plays a key partEil Variations in the lo- 
cal DOS |*KB£f calculated as a function of ferromagnet 
thickness. l3'Ej These-calculations wer£-.all based on the 
quasiclassical UsadelEZl or EilenbergeiEj equations. The 
use of quasiclassical techniques may not be appropriate 
when the thickness of the materials is only of a few atomic 
layers. Also, the Usadel equations are restricted to the 
limit when the mean free path is much smaller than any 
other relevant length scale in the problem, and therefore 
their use in these situations is questionable. It is there- 
fore desirable to study finite-sized systems using a micro- 
scopic, self-consistent theory that can accurately account 
for these geometrical effects. 

We are aware of no work that addresses the influ- 
ence on the proximity effect of the mismatch between 
the three Fermi energies (or Fermi wave vectors) present 
in Fj S junctions (corresponding to the superconductor, 
and to the jn*> and down spin bands in the magnet). Pre- 
vious workEj on this question was limited to the case 
of nonmagnetic metal, at temperatures near T c . It was 
found that when the Fermi wave vector in the normal 
side is smaller than that in the superconductor, a strong 
suppression of the pair amplitude in the normal metal 
ensues. p. Also for nonmagnetic materials, the DOS was 
studiedEI using a microscopic formalism that allowed for 
different Fermi wavevectors, in the context of layered 
short-coherence length superconducting structures, but 
there was no systematic study. Therefore this influence 
is still an open question in the full parameter range. For 
F/S junctions, the influence of Fermi wave vector mis- 
matches on the proximity effecjUas virtually uncharted 
territory. Spectroscopy studiesEJ revealed a nontrivial 
dependence of the conductance spectrum on Fermi wave 
vector mismatch, however the proximity effect was ig- 
nored there and the calculation was not self consistent. 

Another relevant quantity that has a strong influ- 
ence on the proximity effects and which has been in- 
sufficiently studied, is the interfacial scattering strength. 
The variation of T c with interface scattering strength was 
calculated,u and the influen^r©!' interface scattering was 
investigated experimentallyBEia for S/F/S structures. 



High-sensitivity transport measurementsES revealed that 
interface barrier strength was an important parameter. It 
is thus desirable to study the effects of varying degrees of 
barrier strength on the characteristic proximity lengths. 
Since appreciable scattering at the interface should lead 
to a reduction in F(r) near the interface, this is another 
example where a systematic, self-consistent solution to 
this problem is needed. 

The main aim of this paper is therefore to present 
an extensive theoretical investigation of the influence of 
these many relevant parameters on the F j S proximity 
effects. We will use for these purposes a very recently 
developed numerical method that allows for the exact 
self consistent solution of the relevant microscopic equa- 
tions. The approach is based on numerically solving, the 
continuum Bogoliubov de Gennes (BdG) equationsa for 
the quasiparticle amplitudes. The method of numeri- 
cal self consistent solution has been described in Ref. |34| 
where results for particular cases (zero temperature, no 
barrier or mismatch, and semi-infinite geometry) were 
given. These procedures allow for the study of coherence 
lengths much longer thap-,those one can .consider in lat- 
tice real-space models.E3E3. It was shownEJ that for F/ S 
junctions there is, besides the usual characteristic spatial 
period £ 2 , another length scale £1 « £2 which describes 
the fast decay of F(r) very near the interface. The above 
calculations were performed only at T = 0. Finding the 
temperature dependence of the F/S proximity effects will 
thus be a part of our task here. 

Our objective is to investigate the length scales charac- 
terizing the F I S proximity effects in both bulk and finite 
sized junctions consisting of a ferromagnet of varying po- 
larization, (including the nonmagnetic limit) in contact 
with a superconductor. As alluded to above, the often ex- 
treme differences in length scales in the problem require 
a self consistent microscopic theory that can deal with 
them simultaneously without the approximations inher- 
ent to quasiclassical and dirty-limit equations. We shall 
explore the whole of the parameter range including the 
effects of temperature, Fermi energy mismatch, interfa- 
cial scattering, and finite sample size. Results will be 
given for the pair amplitude and for the local DOS for 
both bulk and finite heterostructures, and thus we will 
analyze the various length scales involved. 

Although the objective of this comprehensive study is 
to stimulate new experiments and to help analyze and 
characterize the resulting data, we aim also to make 
contact with existing experimental work. Recent STM 
measurements^ indicate a clear modification to the nor- 
mal metal density of states for a Nb-Au junction, as a 
function of superconductor width. We compare our re- 
sults with this data using a bilayer model in which both 
the normal metal and superconductor have widths of or- 
der £o- We also compare our theoretical results with tun- 
neling dataEl from a F/S (Ni-Al) junction in which local 
DOS measurements were taken in the superconductor. 
The relatively large exchange energy of Ni makes it an 
ideal candidate for investigating the effect of magnetism 
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on the pairing correlations in the superconductor. In 
both cases, we find, using relevant values of the param- 
eters describing the materials used and the geometry of 
the experimental set up, very good agreement between 
theory and experiment. 

The rest of this paper is organized as follows. In Sec. ||, 
we outline the method of self consistent solution to the 
problem. In Sec. Ill, we present our results for the numer- 



ous parameters discussed above for different geometries, 
and compare our results with recent tunnel spectroscopy 
data. Finally in Sec. |fy] we summarize our conclusions. 



II. METHOD 

In this section we briefly outline our basic equations 
and methods. We begin with a brief review of the spin 
dependent microscopic BdG equations in our geometry, 
and then outline the numerical method used for solv- 
ing them. We also explain the procedure for calculating 
physical quantities paramount in the study of proximity 
effects, namely the pair amplitude and the local DOS. 
Most of the techniques used here follow those of Ref. 
|34]. We will omit many of the details given there and 
we focus our attention below on those points where the 
methods employed depart from those developed in that 
work, such as the inclusion of an insulating barrier and 
of finite temperatures. n 

The BdG equationsa are a conceptually simple and 
convenient set of microscopic equations used for study- 
ing inhomogeneous superconducting systems, in our case 
structures involving, in addition to the superconductor, 
a ferromagnet or a non-magnetic normal metal. We con- 
sider a particular slab-like geometry where the materials 
are assumed to extend to infinity in the x — y plane, and 
have a total arbitrary thickness d along the z direction, 
where the only geometrical variation occurs. We denote 
the thicknesses of the ferromagnetic and superconducting 
layers by d! and d—d' respectively. These materials are in 
general separated by a thin insulating barrier at z = d'. 
Since in this geometry the system is translationally in- 
variant in the x — y plane, some aspects of the problem 
are effectively one-dimensional. For this configuration, 
we then have two sets of coupled equations, one for the 
spin-up and spin-down quasiparticle and quasihole wave 
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, and another for (u^v^). The first 

(h = k B = i), 

+ U(z) - E F (z) - h (z)]ul{z) 

+A(z)vi{z) = e n ul(z), (la) 

. + U(z) - E F (z) + &o(*)]«i(z) 

+A(z)ul(z) =e n vi(z), (lb) 



where sj_ is the transverse kinetic energy, e n are the 
quasiparticle energy eigenvalues (the index n labels the 



relevant quantum numbers), ho(z) = hoQ(z — d') is the 
magnetic exchange energy. Scattering at the interface 
(assumed to be spin independent) is accounted for by 
the potential U(z) — H8(z — d'), where H is the bar- 
rier strength parameter. The pair potential A(z) satis- 
fies a self-consistency condition as discussed below, and 
since we will assume that there is no current flowing in 
the system, we can take it to be real. In general, we 
must allow for the possibility of having m to three dif- 
ferent Fermi wavevectors or band widthsE] in the prob- 
lem. The quantity Ep(z) equals EpM in the magnetic 
side, < z < d', so that Ep-\ = Efm + ^0, an d 
Epi = Epm — ho, while in the superconducting side, 
d! < z < d, Ef(z) = Efs- We will be assuming parabolic 
bands so that e± = l/2m(k 2 + k 2 ) and there are three 
Fermi wave vectors, corresponding to Ep^, Ep± and Eps- 
The solutions for the other set of wavefunctions (u^vjj 
are easily obtained from those of Eqns. (|l|) by allowing 
for both positive and negative energies, and then using 
the transformation: — > v\ , — * —u^, e„ — * — e n . 

Equations (|l|) must be solved in conjunction with the 
self consistency condition for the pair potential, 



A(z) 
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E 



K(z)vl(z)+ul(z)vl(z)] tanh(e„/2T), 



(2) 



where T is the temperature, and g{z) is the effective BCS 
coupling constant, which will be taken to be zero outside 
the superconductor and a constant within it. The prime 
on the sum in (^J) indicates that the sum is restricted 
to eigenstates with |e n | < top,, where ujd is the Debye 
energy. 

We now solve Eq. ([!]) by expanding the quasiparticle 
amplitudes in terms of a complete set of functions <f> q {z) 1 
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We will use the normalized particle in a box wavefunc- 
tions, 4> q {z) = y/2/dsm(k q z), as our choice for the com- 
plete set. Here k q — q/ird, and q is a positive integer. 
The finite range of the pairing interaction ujp, permits the 
sums in (^) to be cutoff at an integer N as discussed in 
Ref. |34|, in a way that depends on the maximum wavevec- 
tor present. Upon inserting the expansions (||) into (Q), 
and making use of the orthogonality of the <f> q (z) , we ar- 
rive at the following 2N x 2N matrix eigensystem, 



H+ D 

D H~ 



(4) 



where ^ n is the column vector corresponding to 
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). The matrix elements are 



dz (f> q (z)U(z)(j>q>(z) 
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pd' pd 

E F i I dz <fr q (z)(j>qi (z) — Eps \ dz<j) q (z)(f> q ,(z), (5a) 

JO Jd' 

d 

6„„> - I dZ<f> q (z)U(z)<f>qt(z) 



TABLE I. Dimensionless variables 



H qq , 



2m 



■ e± 



+E Fl dz <fiq(z)<f)q'(z) + Eps / dz(f>q(z)4> q >(z), (5b) 

J Jd' 



D qq , 



dz(j)q{z)A(z)(j>q>(z). (5c) 



The self-consistency condition is now transformed into, 



A(z) = ^Y,Y, tanh( £n /2T) 

p,p' n 



(6) 



where the sum over the quantum numbers n encompasses 
a sum over the continuous transverse energy e± and a 
quantized longitudinal momentum index q, 



E'-E'E' 



(7) 



The matrix eigensystem Eqn. (Q) and the self- 
consistency condition (||) constitute the primary equa- 
tions drawn upon in this paper. They are solved numer- 
ically, using the algorithm developed and described in 
previous workEJ. The iterative computational process is 
completed when the maximum relative error in A {z) be- 
tween successive iterations is less than a prescribed value 
as explained below. 

Once we have the self-consistently calculated eigenval- 
ues and eigenvectors, we can then construct all relevant 
physical quantities. For example, the usual penetration 
depths are conveniently obtained from the pair amplitude 
F{z), 



F(z) = A(z)/g(z). 



(8) 



The pair amplitude, unlike A(z), is therefore not re- 
stricted by the coupling constant to vanish in the non- 
superconductor. F{z) gives a quantitative measure of the 
superconducting correlations in both the superconductor 
and non-superconductor where there may exist phase co- 
herence between particle and hole wave functions. The 
value of F{z) in the non-superconducting region however, 
does not affect the quasiparticle dynamics since it is only 
A(z) that enters into the BdG equations. 

We can also use our numerical results for the excita- 
tion spectra to calculate the experimentally accessible lo- 
cal single particle properties via the thermally broadened 
density of states (DOS) 

N(z,e) = N 1 (z,s) + N l (z,s), (9) 

where the local DOS for each spin state is given by@ 



Physical quantity 



Dimensionless form 



Exchange energy 

Fermi wavevector mismatch 

Temperature 

Coherence length 

Debye energy 

Barrier strength 

Distance relative to interface 



I = ho/ Efm 
A = (k FM /k FS f 
t = T/T c 
Ho EE kFs£,o 
LO EE UJd/EfS 

H b = mH/kFM 
Z ee kFs(z — d) 



7V T (Z,6) -^^["Wp'^^'W^" 6 ") 
p,p' n 

+vi p vl p ,<t> P (z)4> P '(z)f(£ + (10a) 



p,p n 



+ e »)]- ( 10b ) 

Here /'(e) = df /de is the derivative of the Fermi func- 
tion. We will also be interested in the quantity 



5N(z,£) = Ni(z,e)-N l (z,e), 



(11) 



which will be used to characterize the effective leakage of 
magnetism into the superconductor. 



III. RESULTS 

In this Section, we present numerical results for the 
pair amplitude and local DOS, and discuss other phys- 
ically meaningful quantities arising from the self con- 
sistent excitation spectra. We will analyze the various 
length scales characterizing the influence of the super- 
conductor on the ferromagnet and vice versa. Since we 
will consider systems with a wide range of superconduc- 
tor and ferromagnet widths and physical parameters, we 
divide this section into four different subsections dealing 
with the following topics: (I) systematics of the temper- 
ature, exchange field, Fermi wave vector mismatch and 
barrier height for bulk F/S systems. (2) dependence of 
the results on the finite thickness of cither the F or the 
S layer, and finally (3) a comparison with experimental 
results. 

Most of the results are conveniently expressed in terms 
of the dimensionless quantities compactly defined and 
listed in Table |. Unless otherwise indicated, we use 
lu= 0.1 for the Debye cutoff in units of Eps and So = 
kps^Q—50 in this work. All lengths z are measured in 
units of the inverse of kp$. For example the widths of 
the ferromagnet or superconductor layers are represented 
as, Dp = kpsd', and D s = k FS (d — d'). The bulk case 
is studied by choosing values of Ds and Dp sufficiently 
large so that the results become independent of these 
quantities. 
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FIG. 1. The pair amplitude F(z), normalized to the zero 
T bulk value Fo — Ao/g in the superconductor, plotted as a 
function of dimensionless distance Z — kpsz from the inter- 
face. The top panel depicts the normal metal (I — 0) region, 
while the bottom panel shows the superconducting region. In 
both cases the curves correspond, from top to bottom, to tem- 
peratures t = T/T c = 0,0.2,0.4,0.6,0.8,0.9,0.94,0.98,0.99. 
Note the different vertical and horizontal scales used in both 
panels. 



energy cutoff. The value 5000 is five to ten times larger 
than that used in previous work, which makes for bet- 
ter convergence and smoother results. Self consistency is 
achieveda via an iteration process. The process termi- 
nates when the relative error between successive A(z) is 
less than a suitable number, chosen here to be 10 -4 (one 
tenth of the criterion in previous work) . The pair poten- 
tial settles down, after starting with a step function initial 
guess, to its self-consistent form within about twenty five 
iterations. This value is typical for most parameter val- 
ues and system sizes used in this paper. The only excep- 
tions are when the temperature approaches T c (the bulk 
transition temperature of the superconductor), or when 
the superconductor width is of order of Then the 
number of iterations needed for self-consistency is much 
larger (up to several hundred) if one starts with a step 
function guess for A(Z). This problem can be alleviated 
by calculating A(z) self-consistently for a given tempera- 
ture and then use this as input for a nearby temperature 
as described above. 



A. Systematics of the Parameter Dependence 

We consider in this Subsection the dependence of the 
results on temperature and on material parameters (ex- 
change field, wave vector mismatch and barrier height), 
in the limit where both ferromagnet and superconduc- 
tor are very thick ("bulk" limit). By this we mean that 
both Dp and D$ are taken to exceed the temperature 
dependent BCS coherence length. In this subsection we 
have taken Dp — D$ — 12So (recall So = 50), so that 
Dp,D$ 2> So and we are in the bulk limit except ex- 
tremely close to T c , t < 0.99. We subdivide the analysis 
into three categories that address respectively tempera- 
ture and exchange effects, Fermi energy mismatch, and 
interfacial scattering effects. 



As mentioned above, we employ a numerical algorithm 
as in Ref. IIP to solve the self consistent eigenvalue prob- 
lem Eqns. (p), (^|). The procedure involves making a rea- 
sonable initial guess for A(z), where the coordinate z is 
of course discretized for numerical purposes. The initial 
guess may be taken to be a previously obtained converged 
result corresponding to a slightly different set of parame- 
ter values, or, in the absence of any such suitable previous 
result, a step function, e.g., A(Z) = A Q(Z), where Z is 
the dimensionless distance from the interface (see Table 
||), and Ao the bulk value of the gap at T = 0. We then 
diagonalize the 27V x 27V matrix described by Eqn. (||) for 
each value of Ex- The cutoff number TV, as explained in 
Ref. |34| depends chiefly on D = kpxd 7 where kpx is the 
largest Fermi wavevector in the problem and d the to- 
tal thickness, which is up to 1200fc^g in our calculations. 
We use 5000 different values of e± consistent with the 



1. Temperature and exchange energy dependence 

We now consider how variation of the temperature af- 
fects the pair amplitude and local DOS for select values 
of the dimensionless exchange energy /. To isolate these 
effects, we assume that there is no interface barrier and 
no Fermi energy mismatch (A = 1, Hb = 0, see Table |). 
We first examine the case where the non-superconductor 
is a normal metal (1 = 0). In Fig. [l], the pair amplitude 
F(Z) is shown in the normal metal and superconductor 
sides respectively for a wide range of temperatures. In 
all plots, we normalize F(Z) to the zero T bulk value 
F = Ao/g. The two regions Z > (superconductor) 
and Z < (normal) are plotted in separate panels be- 
cause their significant features occur over different length 
and vertical scales. The pair amplitude however, is con- 
tinuous across the interface. We see that in the normal 
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FIG. 2. Normalized pair amplitude as a function of dis- 
tance, in both the ferromagnetic (top panel) and the super- 
conducting side (bottom panel), at the same temperatures as 
in Fig. [lj. The data is for I = 1/4, with all other parameters 
as in Fig. [j]. In the top panel, the amplitude of the oscillations 
decays monotonically with increasing t. 

metal, Fig. [j] (upper panel), F(Z) has a different func- 
tional form at zero temperatures (top curve) than at fi- 
nite temperatures. At T = F(Z) has a very slow decay 
into the normal region, and is expectecB to decay as the 
inverse of the distance from the interface, 

F(Z) = (12) 
\Z\ + c 2 

where ci and C2 are constants. We find that the expres- 
sion Eqn. (|l2| ) is valid, but only in a fairly narrow range 
close to the interface. The actual behavior over the range 
shown is more complex, which reflects that the decay of 
F(Z) takes place, as we shall see, over two length scales. 
Upon increasing T, thermal effects reduce the phase co- 
herence of the electron- hole wavefunctions and. the rele- 
vant length scale in the normal metal is set byQ 

£ N (T) = v FM /2nT. (13) 

where vfm = Ufm/tti. It is clear from the top panel 
of Fig. | that as T increases the length scale character- 
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FIG. 3. Normalized pair amplitude for I — 1/2, plotted as 
in Fig. ^, for the same parameter values and temperatures. 
Note the reduction in the characteristic period compared to 
Fig.| 

izing the decay of F(Z) decreases. The decay of F(Z) 
at a fixed, finite temperature cannot be fit to a single 
exponential in all of the spatial region shown. For tem- 
peratures close to T c , and far from the interface, an ap- 
proximate form for the pair amplitude has been givenp 

F(Z) = $(Z) exp(-\Z\/£ N (T)), (14) 

where <&(Z) is a slowly varying function. Our results 
agree with Eqn. ( |l4| ) in the temperature regime near T c 
(t > 0.9), and for sufficiently large \Z\. (\Z\ > So). 
However, the overall decay of F(Z) is more complicated 
and cannot be described by a single exponential decay. 
There always seems to be a second length scale in the 
problem, even near T c . 

Turning now to the superconductor side, the lower 
panel of Fig. |l| shows the normalized pair potential F(Z) 
at the same values of T used in the panel above it. At 
T = 0, the characteristic decay length of the pair ampli- 
tude is given by the usual zero temperature BCS coher- 
ence length £o- As the temperature is increased however, 
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the depletion of superconducting correlations occur over 
a length scale which increases with T. We denote this 
scale (in units of inverse &fs) by £s(T). For tempera- 
tures close to T c the profile for the pair amplitude is well 
known from standard Ginzburg-Landau theorytj, and has 
the following form 



0.05 



F(Z) = F (T)tanh 



Z + Z a 



\/2&(T) 



(15) £ 



where Zq is a parameter to be determined by the condi- 
tion 



1 dF 
~F~dz 



z=o 



(16) 



and b is an extrapolation length, which in our dimen- 
sionless raits is of order Sq. These expressions hold 
provided^ that S is not too small. Our results fit 
Eq.(|l5|) adequately for temperatures T < T c over the 
entire Z range shown. As the temperature is decreased, 
the expression in Eqn. ([l5|) ceases to be correct for the 
entire range of Z, but remains an adequate fit within 
a region of at least one coherence length from the in- 
terface. We use this expression, therefore as a fitting 
function even at lower temperatures in order to extract 
the dimensionless length scale £s characterizing the de- 
cay away from the interface. We find that for most of the 
temperature range, the characteristic length £s fits well 
to the Ginzburg-Landau expression 



£a(T) = 0.74 S 



l-t 



1/2 



(17) 



After having shown the dependence of F(z) on inter- 
mediate temperatures < T < Tc for 1 = 0, and having 
verified that our limiting results are in agreement with 
previous theory and expectations for the standard (non- 
magnetic) proximity effect, we turn to the more inter- 
esting case where the exchange energy parameter / is 
finite. We found above that when 1 = 0, the supercon- 
ducting correlations extend well within the normal metal 
at T = 0, but decay more rapidly as the temperature is 
increased. When an exchange field is present, the spin 
degeneracy that existed for I = 0, is removed. The re- 
sult is that the Fermi wave vectors of the spin up and 
spin down electrons, fc^f, kpi, are different, and conse- 
quently a Cooper pair entering the ferromagnet acquires 
a net center of mass momentum. The superconducting 
order induced in the ferromagnet arises from the prod- 
uct of particle and hole wave functions (e.g., u} l (z)vj 1 (z)) 
summed over all quantum numbers n. It is the super- 
position of these wavefunctions that, causes the super- 
conducting wavefunction to oscillateO on a length scale 
set by the difference in the spin up and spin down wave 
vectors, £2 ~ (&ft ~~ ^fj,) -1 - We have 



1 



(A(l + /))i/2 



(18) 
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FIG. 4. Normalized pair amplitude for 1=1, plotted as in 
Figs. I and §, for the same parameter values and tempera- 
tures. Note the deviation from the simple sinusoidal decay of 
Eqn. ^TJ on the magnetic side. 



where A iSr-the wavevector mismatch parameter of Ta- 
ble |. Sinceo kp^2 ~ 1/1, we see that the characteristic 
length of oscillations scales as 1/1, and therefore, except 
at extremely small /, it is much smaller than length scale 
set in the normal metal case above. 

We have previously studiedo the explicit form for the 
pair amplitude in the ferromagnet at zero temperature. 
We found that for most exchange fields, and except ex- 
tremely near the the interface, the pair amplitude is given 
by 



F(Z) 



sin[Z/(fc FS 6)] 

* {z/(k FS b)) ' 



T = 0, 



(19) 



where a is a constant. Very close to the interface, the 
pair amplitude monotonically decays over a characteristic 
length ^1 which is defined as the first point inside the 
ferromagnet at which F{Z) is zero. The length scale £1 
goes also as kFs£i ~ 1-1 1- 

We are interested in how the pair amplitude and var- 
ious characteristic lengths associated with it are modi- 
fied as T increases, at finite /. We therefore display in 



7 



Fig. g the pair amplitude at both the ferromagnet and 
superconductor for 7 = 1/4 and the same temperature 
values used in Fig. [l| As in that case, the split panel 
arrangement is required by the difference in vertical and 
horizontal scales, but the function F(Z) is always con- 
tinuous at the interface. We focus first on the ferromag- 
netic region. Starting with the top curve in the upper 
panel of Fig. § (T = 0), we see that, beyond a small 
region of fast decay at the interface, the pair amplitude 
exhibits damped oscillations, with a temperature inde- 
pendent period that coincides with the expected value 
kps& ~ 1/^ = 4, independent of T. The envelope decay 
of the oscillations varies inversely with distance as given 
in Eqn. The quantity £i is also independent of tem- 

perature, since as can be seen in the Figure, the location 
of the first node of F(Z) as it monotonically goes to zero 
near the interface is the same for all temperatures. As 
T increases, however, the amplitude of the oscillations in 
F(Z) markedly decreases. This decrease, as we shall see 
below, is not merely a reflection of the smaller value of 
A(T) in the bulk superconductor. Because of this com- 
petition between thermal and exchange energies, the pair 
amplitude now has a slightly more complicated functional 
form than that given by Eqn. (|l9|). We find that in order 
to fit our numerical results, Eqn. ( |l9| ) must be modified 
by incorporating additional spatial and temperature de- 
pendent factors. The amplitude of the oscillatory decay 
of F(Z) no longer decays as the inverse distance from the 
interface, but now has an additional slowly varying expo- 
nential term $'(Z), and a purely temperature dependent 
amplitude, A(t), 



F{Z) = A(t)&(Z) 



sin [{Z + 0)/k FS £ 2 



(20) 



where 8 is a small, weakly 7 dependent shift that accounts 
for the sharp monotonic decay right at the interface into 
the ferromagnet. We find that Eqn. ( f20|) holds for nearly 
the entire range values of 7 < 7 < 1. Certain exceptions 
occur in the extreme cases of very small I ~ Aq/Epm or 
very large 7 ~ 1 and will be addressed below. The tem- 
perature dependence of the amplitude A(t) in Eqn. (20), 



is fitted well by the form A(t) = A(0)(1 - 1 2 ). Thus A{t) 
decreases faster with temperature than the bulk A(T), 
which shows that the decrease of the amplitude with tem- 
perature is not merely a normalization effect but involves 
an intrinsic decrease of the pairing at the interface. Tem- 
perature has a marked effect on the amplitude, but it 
does not wash out the oscillations themselves, which re- 
main quite well defined even at temperatures quite close 
toT c . 

The superconductor side (bottom panel of Fig. |^) 
shows a behavior of F(Z) very similar to that found in 
the 7 = case, with the variation of F(Z) again oc- 
curring over the length scale £s(T). The effect of the 
the exchange field on F(Z) in the superconductor re- 
gion therefore seems to be minimal at all temperatures. 
We will see below however, that the pair amplitude is 




FIG. 5. Variation of the minimum value of the self consis- 
tent pair potential (Amin) in the superconductor as a func- 
tion of dimensionless temperature for (from top to bottom) 
7 = 0,0.25,0.5,1. 



only partially useful in conveying the total effect of mag- 
netism leakage into the superconductor. The quantity 
SN(z) from Eqn. (11) will be used below for extracting 
additional useful information on this question. 

Continuing with larger values of 7, Fig. || shows the 
case of I = 1/2. We first address the ferromagnet side in 
the top panel. F(Z) follows the form given in Eqn. 
with a damped oscillatory behavior similar to that in the 
7=1/4 case. However, one should note in making the 
comparison that the horizontal axis scale here has been 
reduced by a factor of two with respect to that in Fig. || 
since the characteristic period has approximately halved 
in accordance with Eqn. (|l8|). In the superconductor side 
(Fig. H bottom panel), the decay away from the interface 
is governed by the length £s(T). The very slight wiggles 
in F(Z) which may be observed near the interface are 
due to the increased mismatch of the two Fermi energy 
levels in the ferromagnet with Eps- Overall however, the 
superconducting region shows little change compared to 
the previous case. 

We now turn to the extreme (half metallic) case of 
7=1, where only one spin band is present in the ferro- 
magnet at the Fermi level. We plot results in the same 
way, and for the same temperatures as in Figures ^, |^, 
and |§|. The top panel in Fig. ^ illustrates the pair ampli- 
tude in the ferromagnet. The characteristic length scale 
that describes the main oscillatory behavior is given from 
Eqn. (O) as kps^2 = 1/V2, (recall that we are using 
A = 1 in this subsection). The relevant spatial varia- 
tions occur now only on an atomic scale, (see horizontal 
axis). This reflects that Andreev processes are inhibited 
by the absence of Fermi level down states deep within 
the magnet. We also see clear deviations from the pure 
damped sinusoidal behavior seen for the previous two ex- 
change field values. The superconducting region (bottom 
panel) follows the same pattern as the other cases, but 
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FIG. 6. Local DOS (normalized to its Fermi level value in 
the normal state of the superconductor) versus the dimension- 
less energy e/Ao at / = for t = 0.02 (top panel) and t = 0.1 
(bottom panel). The curves shown, from top to bottom at 
e/A < 1, are for for Z = -100, 100, 200 respectively. 

here we see that very near the interface there exist small 
oscillations of order of the Fermi wavelength. The oscil- 
lations were barely glimpsed at / = 1/2 and disappear 
with decreasing /. 

The pair amplitude at the interface at constant tem- 
perature decreases markedly with /, while at constant / 
it decreases with T. We illustrate this in Fig. |[ where 
Amin is the minimum value of the normalized A(Z) in 
the superconductor. This minimum occurs right at the 
interface, and because of the relatively wide horizontal 
scale in the bottom panels of Figs. |l|,^, and it is not 
really possible to read its value from these Figures. We 
see in this Figure that the effect of the exchange field 
is quite pronounced, and that as the temperature ap- 
proaches T c , all curves tend to collapse into a nearly 
straight line tending to zero. The depletion of supercon- 
ducting correlations with exchange field at the interface 
is also quite evident from the data shown. 

Having studied the spatial dependence of the super- 
conducting correlations, it is now pertinent to examine 
the local density of states at various positions on both 
sides of the interface. The local DOS gives further direct 
insight into the proximity effect, and more important, 
it is an experimentally accessible quantity. Its calcula- 



tion is achieved through Eqn. (g) and the computed self 
consistent spectra. 

We again consider first the case I — 0. The top panel 
in Fig. g shows the normalized (all results for the DOS in 
this work are presented normalized to the superconduc- 
tor's normal state DOS at the Fermi level) local DOS at 
t = 0.02, while the panel below it corresponds to t = 0.1. 
For both cases the numerical results are plotted for the 
local DOS at the positions Z — — 2So,2So,4So that is, 
one position in the normal metal and two in the super- 
conductor. Focusing on the very low temperature case, 
(top panel), we consider first the normal metal region 
at Z = —100 (top curve in the subgap region). We see 
there the sawtooth-like pattern characteristic aLthe de 
Gennes-St. James states as predicted long agoc3. The 
DOS is small but finite at the Fermi energy due to filling 
by thermally excited quasiparticles, and then rises nearly 
linearly at small energies. The Andreev bound states are 
illustrated by the peaks in the DOS. These are due to 
constructive interference of the electron and hole wave- 
functions, as they undergo Andreev reflection at the F/S 
interface and normal reflection at the vacuum-normal 
metal interface at the opposite end of the sample. The 
characteristic energy E c of the peaks is determined by 
adding up the phases for a given trajectoryEil. It can be 
seen that, in agreement with theoretical expectations, the 
first peak occurs at an energy E c ss iTVFM/^d' ~ 0.2Ao, 
while the other peaks occur approximately at multiples of 
2E C . The energy scale E c can be .seen directly in the cal- 
culated self-consistent spectrum.Ea In the energy region 
below the gap (e/A < 1), we find that for nearly lon- 
gitudinal momenta (e± ~ 0), there exists roughly three 
excitation branches, at the same energies as the peaks 
seen in the upper panel of Fig. [| These peaks subse- 
quently broaden due to the numerous quasiparticle states 
with momenta nearly perpendicular to the interface. The 
two other curves in the top panel of Fig. ^| correspond to 
points in the superconductor at Z = 2 So and Z = 4 So. 
The bound states clearly flatten out as one moves fur- 
ther into the superconductor, and there are no longer 
any states at the Fermi energy. The BCS gap becomes 
quite evident at the position Z = 4 So, when only a hint 
of nonvanishing DOS can be seen at energies below Ao. 

Increasing the temperature tends to smear the previ- 
ous lower T results. Figure || (bottom panel) illustrates 
this for the case of t = 0.1, and the same positions as 
the panel above it. In the normal metal, the sharp pro- 
nounced zig-zap pattern is now smoothed into a series of 
humps. The same is true for the DOS in the supercon- 
ductor, where thermal excitations fill regions inside the 
gap, while flattening and spreading out the smoothed 
BCS peaks. 

We now can proceed to study the effects of a finite ex- 
change field / on the local DOS in the superconductor 
and ferromagnet. We begin with the ferromagnet region 
at / = 1/2, and t = 0.02. As already seen in Ref. HJ, the 
effect of / on the local DOS is very drastic. The bound 
state phenomena in the normal side relatively far from 
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FIG. 8. Normalized local DOS at Z — Ho in the supercon- 
ductor and t = 0.02 for various exchange fields as labeled in 
each panel. The left column corresponds to intermediate and 
large exchange fields, and the right column to small values of 
/ as indicated. 



s/A 

FIG. 7. Local DOS (normalized as in Fig. ^) in the fer- 
romagnet for 1 = 1/2 and t — 0.02 (left column) at four 
positions near the interface. The right column corresponds to 
I — 1 and the same temperature and positions used in the 
left column. 

the interface are no longer observable at I = 1/2 be- 
cause the overall decay of superconducting correlations 
takes place over considerably smaller distances. On the 
other hand, this decay takes place now in a nonmono- 
tonic matter, which gives rise to a new set of features 
in the ferromagnet, very near the interface: in the left 
column of Fig. we show the (normalized) DOS at four 
positions at and very near the interface. The influence of 
the oscillatory pair amplitude (Fig. becomes evident 
as we examine the four plots in the this column. The 
subgap structure in the top curve (Z = —3) evolves so 
that maxima and minima become reversed at Z = — 1, 
closer to the interface. Comparing with Fig. ||[ we see 
that the oscillating superconducting order has in effect 
induced oscillations in the local DOS as a function of 
position within the ferromagnet, and that the large ex- 
change field induces noticeable particle-hole asymmetry. 
The length scale at which the DOS flips coincides with 
the characteristic distance kps^ii given in Eqn. |[|). The 
first and third panels in Fig. [7| (left) are separated by 
AZ = 2. Although this Figure depicts results obtained 



for a rather low temperature, the oscillatory behavior is 
never completely washed out by the temperature, as re- 
marked above in conjunction with the discussion of F(Z). 

The same behavior holds true a fortiori in the right col- 
umn of Fig. 0, which displays data for the half metallic 
case I = 1, at the same temperature and locations. One 
important difference between this and the left column 
is the spatial scale at which the DOS oscillations occur. 
Since l/I = 1, the complete DOS inversion should oc- 
cur at points separated by an interval SZ of order unity. 
Indeed, the curves in the right column of Fig. [?] reflect 
this. 

In Fig. U we consider the superconductor side of the 
junction. In the left column we illustrate the local DOS 
one coherence length £o from the interface for three dif- 
ferent exchange fields ranging from I = 1/4 to I = 1, at 
temperature t=0.02. In the top panel (I = 1/4) there 
is a wide U-shape opening for energies |e/A | < 1. The 
opening then starts to get smaller and the curve trends 
upwards for I = 1/2, as seen in in the middle plot of this 
left column. There is also a slight decrease in the num- 
ber of states at the Fermi level (e = 0). These results 
are consistent with those previously obtainedo at zero 
temperature. For the half- metallic case 7 = 1, the left 
bottom panel shows a further slight reduction in states at 
e/Ao = 0, but but there is barely a hint of the asymmetry 
found in the ferromagnetic side. 

It is of interest also to study the case where the ex- 



it) 



CD 
N 



0,25 0,5 0.75 1 1.25 1,5 




0,25 0.5 0.75 1 1.25 1,5 



0,25 0,5 0.75 



1.25 1,5 




0,25 0.5 0.75 



1.25 1,5 




0,25 0,5 0.75 



25 1,5 



e/A 

FIG. 9. SN (normalized to the normal state DOS summed 
over spins) in the superconductor for 1=0.25 (left column) 
and 1=1 (right column) at t = 0.02 and positions (from top 
to bottom) Z = nSo,n =1,4. 



change field parameter / is weak, of the order of Aq/ Eps- 
Possible resonance effects have been predicted to occur 
at these small exchange energiesJl3 Results are shown in 
the right column of Fig. ||. For the values of the pa- 
rameters considered in this subsection, it follows from 
Table | and the BCS relation k F £ = (2E F )/(irA ) that 
Aq/Efs = 0.0127 here. We consider the same tempera- 
ture and location in the superconductor as for the larger 
values of I in the left column of the Figure. We start with 
the case / = Aq/Efm, (top curve in this column). Fo- 
cusing on energies |e/A | < 1, we see a dramatic sharp 
peak in the DOS near the gap edge, and five smaller 
peaks at lower energies. Upon doubling /, (second panel 
from top) the sharp peak structure near e/Ao = 1 van- 
ishes, and there are now four subgap small peaks. In- 
deed, we have found that the sharp peak exists only at 
I = Ao/ Efm- Finally, the bottom curve shows that for 
I = 3Aq/Efm, only three small subgap peaks remain. 
It would be of considerable interest to verify experimen- 
tally the appearance and disappearance of the very sharp 
peak at the gap edge for I = Aq/Efm- 

Examination of the bottom (superconductor side) pan- 
els of Figs, [j], ||, H and [|, shows that the exchange 



field seems to affect the pair amplitude F(Z) within 
the superconductor relatively little over any significant 
length scales. We want to study the possibility, how- 
ever, that the differential local DOS SN(e,z), defined in 
Eqn. ( p"T| ) may show, within the superconductor, mag- 
netic penetration over distances much larger than that 
revealed by A(z). We have a hint that this might be 
the case: the results for A(z) in Ref. [54], exhibited no 
significant dependence on /, while SN(e, z) was appre- 
ciably nonzero within a small region in the supercon- 
ductor near the interface. These previous results were 
obtained for the special case where kFS — fepT> a con- 
dition, which corresponds to an I dependent mismatch 
parameter A — 1/(1 + 7), that may yield results different 
from the case A = 1 considered here. 

We examine in Fig. |^, the normalized 5N(Z) for a field 
parameter values 7 = 1/4 (left column) and 1=1 (right 
column). We use t = 0.02 and choose four locations in 
the superconductor, at Z = nSo,n = 1,4. At the posi- 
tion Z — So, (top panels) there is a clear manifestation 
of the magnetic proximity effect through a nonzero value 
of S N near c/Aq = 1. The effect decreases as Z increases 
and, for I = 1/4 it nearly dies out at Z = 4S , that is, 
after several coherence lengths. At I = 1 (half metal- 
lic case) the effect is more prominent and extends over 
larger distances. However, the integral of SN over ener- 
gies turns out to be always extremely small (as we shall 
see below) at these distances, at which only the self con- 
sistent energy spectral distribution shows magnetic pen- 
etration spin-splitting effects. We see that SN vanishes 
at the Fermi level but the details of this fairly long range 
redistribution of energy states are nontrivial and difficult 
to interpret. Nevertheless, that the effect is larger near 
the gap energy can be readily understood if one recalls 
thato the imaginary part of the wave vector of injected 
quasiparticles below the gap (in a non-self consistent ap- 
proach) vanishes as,-tiie gap edge is approached. 

In previous workt-1 it was found as mentioned above, 
that for 1=1 and no mismatch between the spin up 
and superconductor band (Ep^/Eps = 1), the effect of 
the exchange field on the superconductor was small and 
SN(Z) decayed away over a few atomic distances. How- 
ever the current condition A = 1 implies, at I = 1/4 and 
particularly at I = 1 a considerable mismatch between 
kFs and fepT- We shall study this point in detail below, 
in the context of our discussion of wave vector mismatch 
in general. 



2. Fermi wavevector mismatch 

In the r-test paragraph, we have seen that, (as 
previouslyEJ seen in a different context) mismatch among 
the three Fermi wave vectors involved in the problem (or 
the three band widths) may have a considerable effect 
on the results. As such mismatch is experimentally un- 
avoidable, we now proceed to investigate it in some detail. 
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tions still decay in accordance with Eqn. (|T^), the only 
modification being a mismatch dependent amplitude fac- 
tor g(A), 
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FIG. 10. F(Z) at t = 0.1 and I = 0, plotted vs. Z, for 
values of the mismatch parameter A (see Table |j|) A = 0.1 — 1, 
in increments of 1/10. Top panel is the normal metal region: 
the curves, from top to bottom, correspond to decreasing A. 
The bottom panel is for the superconductor, and curves from 
top to bottom correspond to increasing A. Inset: crossing of 
the curves close to the interface (vertical line) on the normal 
side emphasizing that F(Z) is continuous. 

Thus, we will consider values of the mismatch parameter 
A (Table ||) different from unity. We will still keep the 
interface barrier parameter at Hb = 0. 

In Fig. [To] we show the pair amplitude for 1 = 0, 
t = 0.1, and A varying from the previous case of unity 
down to 0.1 in increments of 1/10. We focus on the situa- 
tion where the bandwidth in the normal metal is smaller 
than that of the superconductor. This is the more com- 
mon situation in F/S structures and in any case, it turns 
out to lead to more prominent effects. The top panel 
in Fig. |l^ shows an overall suppression of superconduct- 
ing correlations with decreasing A (increasing mismatch). 
Extremely near the interface, \Z\ ~ 1, (see inset), F(Z) 
drops rapidly and the curves cross as A increases. At 
A = 0.1 (bottom curve in this main panel) the phase 
coherence is virtually destroyed at the the distance of 
Z = —200. Away from the interface, the pair correla- 



F(Z) 



g(A) 

Z\+C 2 '' 



(21) 



where g(A) is an increasing function of A. Thus, a smaller 
bandwidth in the normal metal tends to restrict the in- 
flux of Cooper pairs. Physically, since the parallel mo- 
mentum of a Cooper pair at interface is conserved, the 
longitudinal component is restricted by the-smaller num- 
ber of states accessible in the normal side.cJ This is con- 
sistent with the bottom panel of Fig. |l^ which shows 
F{Z) in the superconductor. The top curve, correspond- 
ing now to A = 0.1, shows a smaller characteristic length 
of decay from the interface than that for the A = 1 case 
(bottom curve). The subsequent value of A(Z) at the 
interface, Amin, decreases smoothly as A increases. The 
effects of values of A in the range A > 1, (not shown) are 
in the opposite direction, but always much less promi- 
nent. For this reason this range has been deemphasized. 

The general trends are similar in the ferromagnetic 
case. The top panel of Fig. |ll| displays the damped os- 
cillations of F(Z) in the magnet, at I = 1/2. The period 
of the damped oscillations varies inversely with y/A, in 
agreement with Eqn. (|l8|), but they nearly wash out when 
A = 0.1. Also, the sharp monotonic decline very near the 
interface increases in slope with greater A, so that F(Z) 
first reaches zero at a greater distance from Z = 0, thus 
also increasing kpg^i defined earlier. Quantitatively, one 
can obtain an excellent fit for the damped sinusoidal de- 
pendence of the pair amplitude by using Eqn. (|2(]) with 
kFs(,2 as a fitting parameter. The results of doing this 
yield values in excellent agreement with Eqn. (|Tg|). 

In the bottom panel of Fig. [H], the top five curves 
show the drop in F{Z) within the superconductor as the 
interface is approached. The main feature that stands 
out is that the results in the range 0.7 < A < 1 are 
nearly independent of A, while those for A < 0.7 exhibit a 
marked A dependence, similar to that seen in the bottom 
panel of Fig. [h]. This unexpected result arises as at 
7 = 1/2 and A = 2/3 one reaches the special point where 
Ep] = Eps- This property is further exemplified in the 
inset where we present Amin as a function of A. One can 
see a kink in the curve at about A = 2/3. 

The Fermi wavevector mismatch influences also the lo- 
cal DOS as shown in the following figures. Starting again 
with I — 0, we present in Fig. [12] the normalized DOS 
in (left column) the normal metal at Z = — 100 for three 
different values of A ^ 1, and (right column) for the su- 
perconductor at Z = ISo, at the same A values. The 
corresponding A = 1 results are in Fig. g. Both cases 
are at the low temperature of t — 0.02. As A decreases, 
(higher mismatch) we see that, for e/Ao < 1, the bound 
state peaks decrease, until they disappear at A = 0.1. 
The decrease at higher energies reflects our normaliza- 
tion. The superconductor side (right column) shows an 
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FIG. 11. Normalized pair amplitude for the case 7=1/2 
and t — 0.1. The top panel illustrates the variation of F(Z) 
for three different A (as indicated) in the ferromagnet, while 
in the superconductor side (bottom panel) results for the same 
A values used in Figjl^ are shown (in order of decreasing A 
from top to bottom). The inset depicts the value of A(Z) at 
the interface, Amin as a function of the parameter A 



interesting trend. We examine there the point Z = 50, 
one coherence length away from the interface. As one 
decreases A, the small but distinct peaks within the gap 
turn into small wiggles at A = 0.4 and disappear alto- 
gether for A = 0.1. The BCS peaks at the gap edge 
become much more pronounced, indicating a substantial 
reduction in Andreev reflection at the interface because 
of the increased mismatch in Fermi energies, which re- 
sults in superconductivity being more "confined" to the 
superconductor. 



We saw in Fig. 11 that the damped oscillations of the 



pair amplitude inside the magnet increase in wavelength 
and decay quicker with decreasing A (increasing mis- 
match). We investigate the effect this has on the DOS 
in Fig. [[3] for the same value of I, I = 1/2, as in that 
Figure. The left column of Fig. |l^ depicts the changes in 
the magnet side local DOS associated with the same vari- 
ation in A presented for I = in Fig. [L2[ Corresponding 
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FIG. 12. Normalized local DOS plotted versus the dimen- 
sionless energy e/Ao at / — 0, and t = 0.02 for three values of 
the mismatch parameter A labeled in the left panels. The po- 
sition is at Z = —2 So in the normal metal side (left column), 
and at Z = So in the superconductor side (right column). 



results at A = 1 were given in Fig. 0. As expected the 
DOS again experiences oscillations correlated with the 
characteristic length kps£,2, as a function of A. The left 
arrangement of panels in Fig. 13 illustrates this point. 
The coordinate is fixed to Z — —4. One can see an evi- 
dent inversion between the A = 0.8 and A = 1 (Fig. |l3| ) 
cases, whereby the positions of minima and maxima are 
interchanged. The superconductor side is examined in 
the right column, which shows the same A values as in 
the left set of panels, at a distance Z of one correlation 
length inside the material. The top panel exhibits be- 
havior similar to that found for A = 1: the density of 
states within the gap is appreciably nonvanishing, and 
the peaks at e/Ao = 1 are relatively low. The peaks near 
the gap edge for A = 0.4, below are more prominent and 
there is a concomitant decrease in subgap states. This 
is quantitatively different from what we saw at I = 0, 
where there were more subgap states and the BCS peaks 
were significantly sharper. Finally, the bottom panel re- 
veals a near absence of states below the gap, and the 
usual BCS-like peaks at e/Ao = 1- 

that F{Z) in the super- 



It was also seen in Fig. 



11 



conductor decayed away from its bulk value near the in- 
terface in a strongly A-dependent manner. To address 
whether this parameter also affects the spin-splitting in 
the superconductor, we now calculate 8N(z, e) defined in 
Eqn. (|ll|). Figure || shows SN (still at t = 0.02 and nor- 
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FIG. 13. Normalized local DOS for I = 1/2, and with the 
same temperature and A values used in Fig.hi The left col- 
umn corresponds to the position Z — —4 in the ferromagnet, 
while the right column corresponds to Z = Hq in the super- 
conductor. 
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FIG. 14. Normalized SN [see Eq.0, at Z = Ho in the 
superconductor for different values of A. 



positive. As a rule, spin-splitting effects in the self consis- 
tent DOS extend through several times £o (being larger 
near e/A = 1 for the reasons already discussed), except 
of course at the reversal point. The parameter charac- 
terizing the degree of wavevector mismatch is therefore 
important in the study of proximity effects on both sides 
of the F/S interface. 



3. Interface scattering 



malized as in Fig. |9|) as a function of the dimensionless 
energy and at a distance of one £ from the interface, for 
several values of A, at / = 1/2. Starting at A = 1, we see 
an effect reminiscent of what was seen in Fig. ^: there 
is a net negative spin population for e/Ao < 0.85, then 



for larger energies, a greater number of up spin states, 
which decays quickly so that the two spin states equalize 
for e/Ao > 1.5. Next, consider the case where A = 2/3 

(when fcp T = k FS ). In agreement with Ref . |34| the result parameter to unity for this study. In Fig. p|we present 



Up to this point we have considered only transparent 
interfaces. A thin oxide layer at the interface adjoining a 
superconductor and a normal metal or ferromagnet can 
be modeled by a repulsive delta function potential as de- 
fined earlier in Sec. ||. The spin independent scattering 
strength is parameterized in dimensionless units by the 
quantity Hb, defined in Table .|| 

We fix the temperature to t — 0.1 and the mismatch 



is nearly zero for this special value. This value of A is also 
the point at which the F(Z) plots (Fig. 11) start diverg- 
ing with further decreases in A. At this special matching 
point little leakage of magnetism into the superconduc- 
tor occurs. The importance of this crossover point be- 
comes more evident in the remaining curves, where the 
mismatch parameter is decreased to A = 0.4 and then 
to 0.1. The sign of the SN variations with energy is re- 
versed. This pattern, and the relatively large maximum 
and minimum values of SN reflect that the high peaks 
reached by the DOS at these values of A (see Fig. |lj, 
right column), occur at slightly different values for the 
up and down spin bands. Again, the magnetic moment 
at those distances is very small: if one integrates the 
normalized SN over the variable e/Ao the result is of or- 
der 10~ 2 at Z = S , changing sign at A = 2/3. Only 
very near the interface, at values of Z of order unity, 
we find that this integral is larger, and of course always 



the pair amplitude in both the superconductor and the 
normal (/ = 0) metal, in the format of previous Figures. 
We consider five equally spaced values of the dimension- 
less barrier strength ranging from Hb = to a relatively 
strong interfacial scattering barrier Hb = 0.8. First we 
examine the normal metal side, in the top panel. The 
effect of the barrier is quite pronounced, as the pair am- 
plitude still decays slowly into the normal side, but with 
an overall large decrease in amplitude. Also evident are 
Friedel typeO oscillations in F(Z) near the interface: as 
the insulating barrier becomes stronger, the two parts 
of the system become more isolated from each other. 
The pair amplitude in the region shown is adequately 
fit by the functional form of Eqn. ^ but with the pa- 
rameters ci and C2 being both functions of Hb- Within 
the superconductor, the pair correlations within a range 
of order £o from the surface increase with increasing bar- 
rier strength. This is illustrated in the bottom panel of 
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FIG. 15. The normalized pair amplitude at I = 
for different values of the barrier strength (see Table 
Hb = 0,0.2,0.4,0.6,0.8, and t = 0.02. The top panel is 
for the normal metal region. The curves, from top to bot- 
tom, correspond to increasing Hb- The bottom panel depicts 
the superconductor side. The curves from top to bottom are 
in order of decreasing Hb in this region. The inset reveals 
the crossing of the curves near the interface (vertical line) at 
length scales too small to be seen in the main panel. 



Fig. [l5[ where the rise of F(Z) near the interface can 
be seen to sharpen with increasing Hb- The top curve 
(Hb = 0.8) has the least overall variation in the scale 
shown. The oscillations near the interface have the same 
period as in the normal metal, and their amplitude in- 
creases with Hb ■ As remarked in connection with Fig. [l(], 
F(Z) is continuous at the interface, but here however, the 
curves cross in the superconductor very near the inter- 
face. This is illustrated in the inset since this property 
is not visible in the horizontal scale of the main figure, 
where we emphasize longer-range changes. 

The case of a finite exchange field (with / = 1/2) is 
shown in Fig. 06. All other parameter values are the 
same as in Fig. [15. Examining first the magnet side, (top 
panel) reveals that the amplitude of the damped oscilla- 
tions decreases as the scattering potential is increased. 



N 




FIG. 16. Normalized pair amplitude at I = 1/2, for the 
ferromagnet (top panel), and the superconductor (bottom 
panel). All other parameters and curve trends are the same 
as in Fig.|l5|. The inset illustrates the variation of Amin as a 
function of the barrier strength Hb- 



The period is independent of Hb, in agreement with 
Eqn. ([i~8"|). The additional decay of the amplitude of the 
oscillations can be incorporated into Eqn. ( p0| ) through 
a multiplicative factor that decreases linearly with Hb- 
The location of the first node of F(Z) in the magnet is 
nearly unaffected, demonstrating that both the charac- 
teristic length scales &fs£i and Icfs^ are independent 
of interface transparency. On the superconductor side, 
the bottom panel shows how the length scale over which 
F(Z) regains its bulk value from the interface decreases 
as the scattering potential increases. We again see oscil- 
lations in F(Z) for finite values of Hb, near the interface, 
here they are more marked than in the I — case. The 
value of Amin at a given Hb always decreases with i\ 
This is illustrated in the inset at the bottom of Fig. [if], 
where we plot Amin as a function of Hb , for two values 
of /. The difference between the two curves is largest 
at Hb = 0, while they trend closer with increasing bar- 
rier strength. As the barrier becomes very strong the 
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FIG. 17. The effect of a finite barrier on the local DOS: 
Normalized local DOS at I = for finite values of the barrier 
strength H B = 0.2, 0.4, 0.6 and temperature t = 0.02. The 
left column corresponds to Z — —100 in the normal metal, 
and the right column is for Z = Ho- The value of Hb at each 
panel on the right is the same as that in the panel on its left. 

proximity effects become minimal. 

We conclude this subsection with a look at how the 
local DOS is modified due by finite barrier strength. In 
Fig. [I?], we display the normalized local DOS for I = 0, at 
t = 0.02. The left column shows the results for the nor- 
mal metal side. The curves all correspond to Z = —100. 
The Hb = results are in Fig. while here we show 
results for increasing values of Hb as labeled in the fig- 
ure. The general trend on increasing the scattering po- 
tential is a reduction in the magnitude of the peaks for 
subgap energies. The characteristic energy spacing E c 
shows relatively little change, but the shape of the peaks 
is drastically altered. The corresponding change in the 
DOS profile for the superconductor is shown in the right 
column, where we present the local DOS for the same 
values of the barrier strength and at the point Z = So- 
The de Gennes St. James bound states still evident 
in the top curve become smeared out until at the bot- 
tom curve, where Hb — 0.6, the influence of the normal 
metal becomes almost nonexistent. Thus we find that 
although both the insulating barrier and the Fermi en- 
ergy mismatch tend to destroy superconducting order in 
the non-superconductor, their DOS signature is quanti- 
tatively different. 

The effect of barrier strength at finite exchange fields 
can be seen in Fig. ||. We take 7=1/2, with all other 
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FIG. 18. Local normalized DOS for 1 — 1/2 and the same 
barrier strengths and parameter values used in Figjl^. The 
distance on the left side, however, is Z — —4 (left column) 
while on the superconducting side (right column) we have still 
Z = 50. 

parameter values at the same values as in Fig. |l7| . In the 
ferromagnet side, the top left curve in Fig. n8T demon- 
strates a wide structure in the subgap DOS. Upon in- 
creasing the barrier strength (lower panels), there is a 
dramatic reduction in this structure. As can be seen, 
when the barrier strength is rather large (Hb = 0.6), the 
DOS shows very minimal signs of the proximity effect at 
the distance from the interface considered. The super- 
conductor side (right column), for Z = So and the same 
values of Hb as the other panel, reflects the trend seen in 
the ferromagnet. The subgap states at zero barrier (top 
curve, see also Fig. ^) gradually disappear with increasing 
Hb, while simultaneously sharp BCS peaks develop. For 
Hb > 0.5, the results follow closely the 1 = case (see 
Fig. |l7](b)) since the influence of the non-superconductor 
material has vanished. These results illustrate the impor- 
tance of fabricating samples with good, clean interfaces. 

B. Structures 

All of the above results pertained to "bulk" structures, 
in that both slabs were taken to have dimensions signif- 
icantly larger than the zero temperature BCS coherence 
length. We now address what happens when either the 
ferromagnet or the normal metal is thin enough so that 
size effects are appreciable. A bilayer system of this type 
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FIG. 19. Pair amplitude in a superconductor in prox- 
imity to a non-superconducting layer of finite thick- 
ness Df (see text). The top panel shows results for 
1 = 0. The curves from top to bottom correspond to 
D F = 0,5,10,20,30,40,50,100,200, respectively. The bot- 
tom panel has results for I = 1/2. The top two curves are for 
Df = 0, 5. The other curves, all of which essentially coincide, 
are for the remaining values of Df as shown in the top panel. 
The insets show Amin vs Df in each case. 



is an appropriate model for the case when the mean free 
path in the finite layer is larger than the layer's width. 
We will present a broad range of results, varying Dp from 
a few atomic spacings up to of order So, while keeping 
Ds ^S> S , and vice versa. We will con sider the case 

For the 



where both Dp and Ds are small in Sec. [II C 



sake of brevity, we will take the interface to be transpar- 
ent (Hb = 0), the mismatch parameter to be A = 1, and 
fix the temperature to t = 0.02. 

We begin with a normal metal (I — 0) of finite width 
backed by a "bulk" superconductor, taken here to be 
Ds = 16So. The top panel in Fig. |l^ shows the pair am- 
plitude in the superconductor for various normal metal 
widths. We show F(Z) only for the superconductor side 
since the pair amplitude in the finite normal side is cut 
off at different distances. The top curve corresponds to 
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FIG. 20. Local DOS for the same geometry analyzed in 
Fig. |l9| at / = 0. From top to bottom, each set of panels cor- 
responds to Df = 20, 50, 100, 200. The left column shows the 
normal metal (7 = 0) local DOS N(z,e), spatially averaged 
over one Ho as described in the text. The right column is the 
local DOS at Z = 50 in the superconductor. 



a single superconductor slab (zero width for the normal 
metal) , while subsequent curves are for increasing normal 
metal widths ranging up to = 200. The oscillations 
near the interface at zero or very small Dp, are again the 
well-known geometrical Friedel oscillations. The largest 
changes in the pair amplitude near the interface occur 
for < Dp < 20. When Dp > 50, the characteris- 
tic length scale for superconducting depletion is given 
approximately by the coherence length £o- The inset 
displays Amin (as usual, the value of A(Z = 0)) as a 
function of Dp. It is seen that Amin drops rapidly un- 
til about Df = 50 = So and thereafter it decays more 
slowly. 

The same geometry as in the top panel of Fig. |l9| but 
with the finite non-superconducting layer being ferro- 
magnetic (I = 1/2) is considered in the bottom panel 
of this Figure, which again depicts the pair amplitude in 
the superconductor. As the ferromagnet thickness Dp 
begins to increase from zero the pair amplitude drops 
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very rapidly, as in the 1 = limit. A notable distinc- 
tion exists here however. When Dp is larger than about 
Dp = 10, the characteristic length scale over which F(Z) 
rises to its bulk value becomes approximately indepen- 
dent of Dp. This behavior is also seen in more detail in 
the inset of the same Figure where we plot Amin as a 
function of Dp. The decay of Amin occurs nearly en- 
tirely in the region Dp < 10, while for / = it takes 
place over a much more extended range. This of course 
reflects that the superconducting penetration (at low T) 
into the normal metal is very large, while for a magnet 
with /= 1/2 it is characterized by a length of order £2- 
Once Dp reaches that limit, further increases are ineffec- 
tive. 

We have also calculated the DOS for the geometries 
used in Fig. Since the non-superconductor layer is in 
some cases quite thin (Dp < So), the local DOS in the 
normal region exhibits strong oscillations as a function 
of Z. For this reason, we present results for the spa- 
tial average of N(Z, e) over a distance in the Z direction 
equal to the layer thickness Dp if Dp < S , or over one 
dimensionless coherence length, S , if Dp > S . In the 
latter case this average is centered at Z = —Dp + Ho/2. 
We present in the left column of Fig. the averaged 
DOS within the normal (/ = 0) metal for four different 
thicknesses, at t = 0.02. The top left panel corresponds 
to a thin film with Dp = 20 (recall the superconductor 
is in the bulk limit). A clear "mini gap" structure is 
present. As Dp is increased to Dp = 50, a much smaller 
gap remains, and multiple ripples rise to two larger bound 
state peaks. If Dp is doubled again to Dp — 100, the gap 
disappears. Upon increasing Dp further to Dp = 200, 
we see another peak emerge and form the initial stages 
of the sawtooth-like profile seen earlier in bulk systems 
(see Fig. ||). Thus we find that there exists a maximum 
thickness for the normal metal Dp w So, such that, if ex- 
ceeded, the gap in the normal side DOS disappears. The 
observed filling in of the states originates from quasi- 
particles with relatively large momenta parallel to the 
interface (fcj_ « kps). 

In the right column of Fig. ^0|, the local DOS in the 
superconductor at a distance Z = So from the inter- 
face is shown. Here we do not spatially average the local 
DOS, since we are in the bulk regime and the DOS varies 
smoothly. We present the local DOS at Z = So, while all 
other parameters take the values used previously in the 
left column. The top curve (Dp = 20) shows a widen- 
ing of the gap, while the main peaks still remain below 
e/A = 1. The panel below demonstrates the bound 
state peaks being pushed further towards the Fermi level. 
For Dp = 100, the single pair of peaks has moved inward 
even further as a marked through develops at e/Ao ~ 1. 
The effect is more pronounced in the bottom curve, where 
Dp = 4S , and the de Gennes St. James states have 
become smaller than the main peaks that have formed 
near the gap edge, which eventually develop into BCS- 
like peaks deeper within the superconductor. 

We consider next the DOS when the non supercon- 
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FIG. 21. Local DOS as in the previous Figure, but for 
/ = 1/2. From top to bottom each pair of panels corresponds 
to D F = 20,30,40, 50. 



ductor layer in the system is a ferromagnet. We present 
results at t = 0.02, and I = 1/2, as in the bottom panel 
of Fig. [l9|. In the magnetic side, we spatially average the 
DOS over its width Dp as described above. The result 
is shown in the left column of Fig. ^l|. The top curve ex- 
hibits two slightly asymmetric peaks at e/Ao « 0.2. The 
structure seen there is washed out at larger Dp. Thcxti 
is no gap in the ferromagnet DOS shown but we founded 
a mini-gap when Dp is small (Dp sa kps&)- The lo- 
cal DOS at the point Z = So inside the superconductor 
is illustrated in the right column of Fig. [2l]. In the top 
curve (Dp = 20), it is seen that the highermost peaks 
are shifted slightly towards lower energies (e/Ao ~ 0.9) 
compared with the bulk BCS result. At even lower en- 
ergies there is a relatively high number of subgap states. 
Upon increasing the ferromagnet 's width to Dp = 30, 
the coarse structure seen previously becomes somewhat 
smoothed out. On the curve below (Dp = 40), numer- 
ous peaks have returned within the gap, and then dimin- 
ish again in the the last curve for Dp = 50. We find 
therefore that the presence of the magnet next to the su- 
perconductor results in more prominent features in the 
DOS, at smaller Dp values. 

We now reverse the role of the two materials in the 
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FIG. 22. Pair amplitude for a structure consisting of 
a superconductor of finite thickness Ds adjoining a thick 
non-superconductor. The main plot in the top panel shows 
the decay of the pair amplitude in the normal metal (J = 0) 
for values of D s = 200, 100, 75, 50, 45. The bottom panel 
shows the pair amplitude at I = 1/2 and the same geometry. 
The values of Ds are the same as in the top panel, and the 
amplitude of the oscillations decays with decreasing Ds- The 
insets illustrate the behavior of Amin vs Ds in each case. 



bilayer, that is, we consider a very thick "bulk" ferro- 
magnet (we take Dp — 16So as was done above for Ds), 
in contact with a finite superconductor layer. Temper- 
ature and other parameters are as in the previous case. 
In order to study fully the geometrical effects associated 
with varying the superconductor thickness, we shall con- 
sider a wide range of widths Dg, taking S close to the 
lower bound, since the superconductor ceases to maintain 
pairing correlations when Ds ^ So. 

The top panel of Fig. [2^ shows the modification of 
the pair correlations in a bulk normal metal (/ = 0) 
that occur as the width of the superconductor varies. 
The top curve (Ds — 4 So) differs relatively little from 
the situation where both the normal metal and super- 
conductor were in the bulk (compare with Fig. [j]). A 
decreasing trend is followed as Ds decreases. The slow 



decay of F(Z) away from the interface is adequately fit 
by Eqn. (|20| ) for Ds > 1.5 So, the only modification 
being an overall Ds-dependent factor that reduces the 
amplitude. The bottom two curves, corresponding to 
Ds = 0.9 So, So have an even slower decay. The inset 
depicts the corresponding change in the pair potential 
at the interface, Aminj as a function of Ds- This inset 
emphasizes the fast rise in the pairing correlations at the 
interface when Ds is on the scale of So, and it includes 
additional values of Ds not presented in the main figure. 

(bottom panel) we show the damped oscil- 



In Fig. 22 



lations of F[Z) within the ferromagnet (/ — 1/2) for the 
same values of Ds as in the top panel. The main effect 
of changing Ds is to reduce the amplitude of the oscilla- 
tions while their period remains, as expected, the same. 
Their amplitude however, drops very markedly when Ds 
approaches So. This is illustrated in the inset, where we 
display Amin versus Ds- The essential behavior is simi- 
lar to that in the 1 = case in the other panel, whereby 
Amin changes the most for Ds < 1.5 Sq. The overall 
magnitude is reduced, however, by the finite value of the 
exchange energy. 

The local DOS is also sensitive to the spatial extent of 
the superconductor. The left column in Fig. ^3| shows the 
normalized local DOS at Z = —100 in the 1 = limit, 
for several values of Ds- Spatial averaging in this case is 
unnecessary. For Ds = 200 (top curve), it is evident that 
this width is sufficient for Andreev reflection to become 
well established, and hence for the complete formation 
of the structure seen in the bulk case (see Fig. ^). At 
Ds = 2So, even though Ds is still larger than So, the 
DOS profile has changed its shape so that the peaks slant 
in the opposite direction. The magnitude of the peaks has 
decreased slightly overall, but their characteristic energy 
spacing E c remains the same with the exception that the 
largest peak has slightly shifted towards smaller energies. 
The next curve (Ds = 1.5So) shows how this trend ex- 
tends further. Finally, when Ds = So, only a slight hint 
of structure remains. These results reflect that the de 
Gennes St. James peaks arise mainly from Andreev re- 
flection at the normal metal superconductor interface, so 
that when Ds decreases, so does the minimum gap (see 
Fig. [2^ inset in the top panel). Next, we examine the 
DOS in the superconducting side (right column), using 
the same parameters as for the left side panels. We per- 
form a spatial average over one coherence length centered 
around Z = Ds — So/2 (in analogy with the thin mag- 
net case). Beginning with the top curve and proceeding 
downwards, we see a rapid filling in of subgap states. As 
Ds is decreased (lower curves), any remnant of subgap 
states becomes smeared out due to the greater influence 
of the normal metal on the superconductor for smaller 
D s . 

We conclude this subsection by considering a finite 
exchange energy of I = 1/2 in the ferromagnet, with 
all other material parameters being identical to those in 
Fig. 23. As mentioned above, for a bulk superconductor 



juxtaposed to a bulk ferromagnet, the oscillatory pair 
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FIG. 23. Left column: localized DOS at Z = -100 for the 
structure considered in Fig. |2^, at I = and superconductor 
widths as indicated. Right column: local DOS for the same 
structure, averaged over one Ho (see text) from the end of the 
superconductor, for the same D$ as in the left column. 
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FIG. 24. Local DOS for the structure considered in Fig. ^ 
with / = 1/2. Left column: localized DOS at Z — —3 for 
1 — 1/2 and superconductor widths as labeled. Right col- 
umn: local DOS averaged over one So from the end of the 
superconductor for the same values of Ds as in the corre- 
sponding left column panels. 



amplitude in the ferromagnet induces oscillations in the 
DOS as a function of position within the ferromagnet (see 
discussion of Fig. ^) . Here, we wish to examine any modi- 
fications on this behavior that may arise from a finite Dg. 
To this purpose, we present in Fig. |24| (left column) the 
local DOS at the position Z = —3 for several values of 
Dg. The top curve, corresponding to Ds = 4So, shows a 
DOS profile with two rounded peaks near the gap edge, 
while the minimum is at e/Ao = 0. This is fairly sim- 
ilar to what was seen at the same distance in Fig. |?]. 
To understand the behavior of the DOS for smaller Ds, 
we recall the spatial dependence of the pair correlations 
in Fig. There it was found that the oscillations in 
F(Z) did not undergo a change in period as the super- 
conductor width decreased. Rather, there was a smooth 
reduction in amplitude as Ds decreased. This suggests 
that changes in the DOS with decreasing Ds would be- 
have similarly. The other curves in the left column of 
Fig. pil agree with this reasoning: the effect of reduc- 



ing Ds is to lower peaks and raise minima, so that the 
leakage of superconducting order is effectively eliminated 
when Ds = So. 

The panels in the right column of Fig. |24| show the cor- 
responding DOS in the superconductor, averaged over 
a distance of one So in the usual way. The top panel 
(Ds — 4So, so that the average is taken centered at 
Z = 175) shows a clear but broad peak at e/A « 1, 
spreading over a significant energy range. This peak de- 
cays rapidly while shifting to smaller energies as Ds de- 
creases. As Ds is reduced to Ds = 2So, the states with 
energies e/Ao ^ 1 fill in rapidly while the primary peak is 
flattened. This trend continues until, at Ds — So, there 
is hardly any evidence of the previous superconducting 
structure, indicating that superconductivity is nearly de- 
stroyed as the thickness is down to one correlation length. 
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C. Comparison with experiment 

We have seen above that the calculated self-consistent 
results for the physical quantities depend in a system- 
atic way on a number of parameters, some of which are 
related to the materials employed, while others are ex- 
perimentally adjustable. While testing these systematic 
dependences is the task of future experimental work, we 
will, in this subsection, compare already existing exper- 
imental data with our theory. _ffie will use data from 
direct local DOS measurementsEffl rather than results 
for indirectly derived quantities. Although we emphasize 
in this work the case where the non-superconductor is 
a ferromagnet, we will examine also the proximity effect 
when the normal metal is a ferromagnet. We will thus 
compare our calculations with the experimental data of 
Rcf. [3^ where DOS measurements were made from the 
normal metal side, and with data from Ref. where 
local DOS measurements in the superconductor side of 
a magnet-superconductor structure were taken. In this 
way, we test our theory against spectroscopy data ob- 
tained by probing either side of the interface, in the two 
cases where the superconductor is in proximity to either 
a normal metal or a strong magnet. In making our com- 
parisons, it is important to make pertinent choices for the 
applicable input parameters, as will be discussed below. 

Consider first the STM data of Ref. |37], where the 
superconductor (Nb) is in contact with a non-magnetic 
metal ( Au) . The experimental configuration consisted of 
a thin layer of Au of thickness 200A, that capped off a 
superconducting Nb dot. The Nb had a smooth relief 
resulting in a thickness (0 — 400A) that decreased away 
from its center. We model this structure, as done in the 
experimental analysisEH, as a bilayer system comprised 
of a normal metal of constant width on a superconduc- 
tor of varying thickness in the z-direction, in a manner 
similar to Sec. Ill B|. We take kpM = 1A , and we as- 



sume the normal metal to have a width of 200 A. The 
transition temperature is T c = 9 K. We assume the tem- 
perature to be T = 270 mK, a value which is slightly 
higher than the experimental value T = 60 mK. This 
is intended, in the usual way, to account for additional 
smearing effects associated with the finite energy reso- 
lution of the apparatus. The interfacial barrier strength 
parameter (Hb) is taken to be zero, which is appropriate 
for the clean, highly transparent interface used, and the 
Fermi wavevectors in the two materials are assumed to 
be equal. Other parameter values used in our calcula- 
tions are the bulk Nb gap value-.Ao = 2 meV, which is 
close to the experimental valueEH, and the Debye cutoff 
parameter lj = 0.03, the value of this parameter having 
little effect on the results. The modeling of the supercon- 
ductor width is less trivial since the Nb dot in fact varies 
in size not only in the z direction, but also in the trans- 
verse direction. With this in mind, we assume a super- 
conductor thickness varying from 50 to 150 Angstroms. 
Because of possible nonuniformities in the composition of 



the Nb dot, this quantity should be viewed as an effective 
thickness that accounts for any geometrical discrepancies 
between our model and the experimental configuration, 
and which may be interpreted to some extent as a fitting 
parameter. The final physical parameter needed is the 
coherence length £o- This parameter must be identified 
here as an effective correlation length to absorb the in- 
herent effects of disorder in the system. Measurements 
were taken at several points: some, (which were labeled 
as points a— din Figure la of Ref. £37|) were on the flatter 
part near the center of the dot. The others, labeled e — j 
in that Figure, were in the sloping part near the edge. 
Disorder effects are likely to be more prevalent in the re- 
gion in the latter points, where oxidation of the sample 
surface has a more pronounced effect on the supercon- 
ducting order. Therefore we have set at £0 = 100A in 
the region corresponding to the points e — j, while for in 
the region between the points a — d we take £ = 200A. 
These two sets of points were recognized as behaving dif- 
ferently in the original experimental analysisEH. 

The geometr y stu died is not quite that in our earlier 
results of Sec. [II B, where we varied the width of the 



superconductor in contact with an infinite normal metal. 
Here, both the normal metal and superconductor are ef- 
fectively thin, and additionally, the assumed value of the 
coherence length is larger. Therefore, separate compu- 
tations were required. We present in Fig. ^5] the com- 
parison of our results (solid line) to the the experimental 
dataED The DOS is scaled to its normal metal limit, with 
the curves shifted by a constant for illustrative purposes. 
The energy is in the same voltage units as in the ex- 
periment. Inspecting the spectra corresponding to the 
points a — d (labeled in our Figure as they were in the 
experimental work) , Fig. ^5] demonstrates the excellent 
agreement between our results and the data. In the top 
panel, a BCS-like gap is most evident for the position 
the location labeled a, and as the effective superconduc- 
tor width is decreased from 150 Ain a to 120 in d, the 
gap becomes smaller, while the BCS-peaks decrease in 
magnitude. The location of the peaks in the fits and in 
the experimental data are seen to essentially agree. A 
similar trend is seen for the remaining probe locations 
in the bottom panel of Fig. ^5], where the effective co- 
herence length and Dg are smaller. The peaks move 
inwards while the previously empty gap starts to fill in, 
with an approximately linear rise near the Fermi level. 
This DOS profile is observably different from that in the 
panel above, where the subgap DOS had a U-shape com- 
pared to the V-shape here. It is remarkable that the level 
of agreement between theory and experiment is so high, 
in that the location of the peaks as well as the origin 
of the minima in the DOS match well over the entire 
spatial range. Thus we find that modeling this particu- 
lar experimental structure as a bilayer is successful over 
the entire spatial range. In Ref. ^ where it was found 
that a fit to all the data using the Usadel equations was 
not possible and that very different physical assumptions 
had to be used for the U- and V-shaped portions. This 
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FIG. 25. Comparison of STM data from Ref. with the- 
ory. Solid curves: theoretical results. Dotted curves: data. 
Vertical scales shifted for clarity. The labels a through j cor- 
respond to different probe positions. See text for details. 



is unnecessary within the exact theory. 

When the normal metal is a ferromagnet, experimen- 
tal studies on the proximity effect are more sparse. The 
continuing advancement in nanofabrication techniques 
however, has made probing the electronic structure of 
F/S nanostructures experimentally accessible and some 
recent good quality data is available. We compare our 
theory with the tunnel spectroscopy data of Ref. || ob- 
tained through probing the local DOS in an Al super- 
conductor adjacent to a Ni ferromagnet. Modifications 
to the DOS in the superconductor are another important 
aspect of the proximity effect which provides useful in- 
formation regarding the influence of the ferromag net on 
superconducting correlations. As predicted in Sec. Ill A 1 



the local DOS near the interface in the superconductor 
should be substantially modified from the bulk BCS re- 
sult. It is then of great interest to see how our results 
compare with the appropriate experiment. 

To test our theory against the Al-Ni experimental data, 
we must choose a set of parameters appropriate for the 
given bilayer. Nickel itself is not a simple Stoner magnet 
with parabolic bands. A nearly free electron monovalent 
metal having the same saturation magnetizationc3 as Ni, 
is easily seen to be to have a value of I of about 0.5 
and this is the value we will use. We assume, a trans- 
parent interface (Hb — 0) in. accordance with the clean 
interface in the experiment E We also take T c = 1.2 K, 
kFM = 0.5A , and in order to limit as much as possi- 
ble the number of input parameters, we keep the Fermi 
wave vectors the same, A = 1. For thick superconduct- 
ing layers the relevant length which governs the dcple-. 
tion of superconducting correlations near the interfaceo 
scales with £0 whenever £0 1. In the experimental 
work, distances were already given in units of the corre- 
lation length, and this makes it particularly convenient 
to compare with theory. Experimental data is given at 
two distances: one far from the interface and the other 
near to it. The precise distance from the tunneling probe 
to the interface in the second case was somewhat uncer- 
tain, however, mainly because of the finite width of the 
probe. We take this position to be 2£ , which is similar 
to the value estimated in the experimental analysis.B 




FIG. 26. Comparison of experimental DOS data for Ni-Al 
structure™ with theory. The symbols represent data taken far 
from the interface (open symbols) and near it (closed sym- 
bols). The curves are the theoretical fits obtained as explained 
in the text. 



Figure 26 shows a fit of our results (solid lines) to the 
experimental data of Ref. ||. The vertical axis is the DOS 
scaled to the normal state value, while the energies are 
in dimensionless units of e/Ao. We account for the in- 
fluence which single-charging effects have on the DOS by 
convolving the DOS calculated from (ph with the accep- 
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tance function P{E) (as described in Ref. ||) that gives, 
the probability for the junction to absorb an energy Ex3 
The curve with the more prominent peak corresponds 
to the DOS in the bulk, away from the influence of the 
magnet, while the other curve is the DOS two coherence 
lengths from the interface. The open and closed symbols 
are the experimental results, the former ones being for 
the bulk DOS. The procedure employed to obtain these 
fits was the following: first we determined the effective 
temperature (to account for instrument resolution) by 
fitting our results to the experimental DOS near the in- 
terface, (closed symbols, solid curve). This resulted in a 
slightly raised effective temperature T — 980 mK (in con- 
trast with the experimental value of T = 100 mK). Then, 
without any further changes, we calculated the DOS at a 
distance of 4£o from the interface (dashed curve). This 
position should represent well the bulk characteristics in 
the DOS, since as we have seen previously, at that point, 
the influence of the ferromanget on the superconducting 
DOS is minimal. No additional parameters were used to 
obtain this second fit. The results are clearly excellent: 
they have the correct peak positions and relative magni- 
tudes. We therefore find again, good overall agreement 
with experiment in this more difficult case. 

IV. CONCLUSIONS 

We have in this work presented extensive results for the 
pair amplitude and local DOS in heterostructures involv- 
ing superconductors and magnetic materials. These re- 
sults were obtained from numerical, self-consistent solu- 
tion of the Bogoliubov-deGennes equations, without ap- 
proximations. We also discussed the length scales char- 
acterizing the influence of the superconductor on the fer- 
romagnet and vice versa. 

For heterostructures with geometric dimensions larger 
than the relevant intrinsic lengths, we have shown in de- 
tail how variation of parameters such as temperature T, 
Fermi wave vector mismatch A, interfacial and scatter- 
ing strength Hb , (see Table ||) affected the pair amplitude 
and local DOS for a wide range of exchange energies /. 
For 1 = 0, and low T, we find that the pair amplitude 
in the normal metal decays approximately as the inverse 
of the distance from the interface, with the overall pref- 
actor depending on A and Hb- At higher temperatures, 
and for A = 1, Hb = 0, the decay markedly increases, 
and is set by a smaller length scale £n(T). The exact 
spatial decay of F(r) was found to be more complicated 
than a single exponential. On the superconductor side, 
we tested our results for temperatures near T c and found 
agreement with standard Ginzburg Landau theory. We 
also extracted the characteristic length of depletion £s(T) 
for intermediate and lower temperatures, something not 
previously done in a systematic way. The length scale 
(at low temperatures) characterizing the decline of F(r) 
near the interface was found to decrease with smaller A 



or larger Hb- The local DOS correlates with these re- 
sults and displays a functional dependence on T, A, and 
Hb as well. We also systematically investigated the ge- 
ometrical effects associated with the finite size of either 
the normal metal or superconductor and discovered that 
a small gap develops in the normal metal DOS when its 
width is small. This "mini-gap" decreases to zero when 
the metal width is close to the coherence length, i.e., 
D F « S . 

At finite values of /, there are two characteristic prox- 
imity length scales £2 and £1 in the ferromagnet, gov- 
erning respectively the spatial period of the damped os- 
cillations in F(r), and the sharp decay at the interface. 
Both lengths vary approximately inversely with J, inde- 
pendent of temperature. The amplitude of the damped 
oscillations however, decreases with increasing tempera- 
ture. Mismatch of the Fermi energies increases the decay 
of the oscillations. A finite barrier strength has no effect 
on the period either, but does reduce the amplitude by a 
factor that scales linearly with Hb- 

The damped oscillatory behavior in the ferromagnet 
induces a corresponding spatial modulation in the local 
DOS. For a transparent interface, the DOS in the su- 
perconductor side (for / > 1/4), exhibited a reduction 
in the usual BCS peaks, with a finite number of states 
within the gap, the number depending on the exchange 
field and location within the superconductor. For small 
exchange fields of order of / ~ Ao /Eps, a significant 
subgap structure emerged and at exactly I = Ag/Eps, a 
resonance phenomenon occurred in which the BCS peaks 
became significantly enhanced. We found there exists a 
long range spin splitting in the superconductor, extend- 
ing over several coherence lengths. A nontrivial behavior 
of SN was found as a function of I and the mismatch 
parameter A: when the point Ep-\ = Ep$ is crossed, the 
spin splitting becomes very short ranged, as found in Ref. 

Finally, we have compared our results with two sets of 
experimental data for the local DOS, corresponding to 
two different values of / and to measurements taken ei- 
ther from the superconductor or the non-superconductor 
side of the heterostructure. In both cases we found, using 
reasonable values for the material and geometric param- 
eters, very good agreement between theory and experi- 
ment. 

The number of relevant parameters involved is so large, 
and the variety of behaviors so rich, that even an exten- 
sive study such as this one must concentrate on the high- 
lights and leave most of parameter space unexplored. It 
is clear however that the machinery developed here can 
be readily applied to most actual experimental situations. 
We hope that this paper will stimulate future experimen- 
tal work and facilitate the analysis of the resulting data. 
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